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A phase from an adiabatic exchange of Majorana bound states (MBS) reveals their exotic anyonic 
nature. For detecting this exchange phase, we propose an experimental setup consisting of a Corbino 
geometry Josephson junction on the surface of a topological insulator, in which two MBS at zero 
energy can be created and rotated. We hnd that if a metallic tip is weakly coupled to a point on the 
junction, the time-averaged differential conductance of the tip-Majorana coupling shows peaks at the 
tip voltages eV = ±{a — 2nl)h/Tj, where a = 7 r /2 is the exchange phase of the two circulating MBS, 

Tj is the half rotation time of MBS, and I an integer. This result constitutes a clear experimental 
signature of Majorana fermion exchange. 

PACS numbers: 71.10.Pm, 03.65.Vf, 74.45.+C, 73.23.-b 


Majorana fermions are charge-neutral quasiparticles 
obeying non-Abelian exchange statistics that may be use¬ 
ful for quantum computation HHT]. Over the past years, 
they have been predicted to appear in certain condensed 
matter systems, like topological insulators [ai , semicon¬ 
ductors with spin-orbit interations mm, and magnetic 
atom chains mHia, in proximity to s-wave supercon¬ 
ductors, and the search for experimental signatures of 
their unusual features has been intensihed. For exam¬ 
ple, the tunneling conductance between a metallic lead 
and a Majorana fermion shows a zero-bias peak [HUTS] , 
and signatures for this prediction have been identihed 
in experiments [MIS]. Other schemes, such as inter¬ 
ferometry in a Dirac-Majorana converter nnun] and a 
Majorana-mediated Josephson effect [52^25] . have also 
been explored. So far, however, most of the signatures 
that are predicted or observed are attributed to Majorana 
features associated with charge neutrality and fermion 
parity anomaly, and less attention has been paid to find 
signatures from the exchange statistics, although several 
schemes to realize exchange or braiding operations have 
been suggested jS] . 

Exchange statistics of Majorana fermions, which differs 
from that of fermions or bosons, can provide alternative 
routes of detection. An adiabatic exchange of two Ma¬ 
jorana fermions 71 and 72 leads to the transformation 
7i 72 and 72 —>■ — 71 , resulting in exotic phase fac¬ 
tors |7]. In Ref. [53], a long circular topological Joseph¬ 
son junction is considered where the energy spectrum of 
the junction is found to depend on the exchange phase 
and fermion parity. Finding further experimental signa¬ 
tures of the exchange phase would be an important task 
for Majorana fermion detection. 

In this work, we propose a transport experiment where 
the exchange phase of mobile Majorana bound states 
(MBS) can be identified. The setup consists of a Corbino 
geometry Josephson junction on the surface of a topo¬ 
logical insulator. By solving the Bogoliubov-de Gennes 
equation, we show that two spatially separated MBS 
at zero excitation energy appear in the junction-lacking 



FIG. 1: Schematic of an experimental setup for performing 
and detecting Majorana exchange. Thin-film superconductors 
labeled 5i and S 2 are deposited on the surface of a topological 
insulator, forming a Corbino-geometry topological Josephson 
junction. In the presence of two flux quanta [2'l>o with $0 = 
h/{2e)], two zero-energy MBS (red balls denoted by 71 and 
72 ) appear at opposite sides of the junction. Their positions 
can be moved slowly while maintaining their relative distance 
by applying a small voltage across the junction, enabling us to 
perform an adiabatic exchange. The resulting exchange phase 
can be identified in the time-averaged differential conductance 
of a metallic tip weakly tunnel coupled to the junction. 


the hybridization between them-when two flux quanta 
are introduced in the junction. Their positions can be 
moved along the circle by changing the superconducting 
phase difference across the junction, allowing for their 
exchange. If a metallic tip is weakly coupled locally to 
the junction, and the two MBS are rotating adiabatically, 
the electron tunneling between the tip and MBS occurs 
every half rotation period. In this weak coupling limit, 
we find that the time-averaged differential conductance 
shows peaks at the tip voltage eV = ±(a — 2'Kl)h/Tj 
with a = 7r/2 the exchange phase [M]. This provides 
direct experimental signatures of the nature of Majorana 
fermion exchange statistics via a dc charge measurement. 

Rotating MBS .— We first find MBS in a Corbino ge¬ 
ometry Josephson junction deposited on the surface {x-y 
plane) of a topological insulator. The setup is illustrated 
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in Fig. The Josephson junction is formed by thin films 
of inner (S'!) and outer (S' 2 ) s-wave superconductors and 
contains two magnetic flux quanta. For simplicity, we 
assume that the distance between Si and S 2 is zero. We 
consider the case where the radius of denoted by R 
satisfies ^ R A so that the vector potential of the 
magnetic field can be neglected, where ^ is the supercon¬ 
ducting coherence length, and A is the Pearl length [55] : 
neglecting the vector potential is justified provided that 
the flux through the area where the bound states are lo¬ 
calized is smaller than the flux quantum |55| . We also as¬ 
sume a short Josephson junction, i.e., 2ttR ^ Aj where 
Xj is the Josephson penetration depth [571 155] . Then 
the Bogoliubov-de Gennes (BdG) equation with Nambu 
spinors 'P = of excitation energy E 

has the form [111911^ 


R-BdO 


nBdG'i'{r,e)=E'^ir,9), 
f H-fi A{r,e)l\ 
-{A*ir,0)l fi-H )■ 


( 1 ) 

( 2 ) 


Here E[ = ■ P is the single-particle Hamiltonian de¬ 

scribing the surface of the topological insulator, a = 
[o'x, o'y) are Pauli spin matrices, I is the 2x2 unit matrix, 
and p is the chemical potential. The induced gap A(r, 9) 
in the presence of the two flux quanta can be described 
by the position-dependent form as 


A(r, 9) 


Aoe*-^! 0 < r < i?, 

Aoe-* 2»+*<^2 r > i?. 


(3) 


where (j)i and (j )2 are spatially uniform phases in each 
region, and the polar-angle-dependent term —29 is due 
to the flux quanta [55] . 

We solve Eq. Q for the case of /i = 0 to obtain Ma- 
jorana zero-energy states, 5 'mi and which are spa¬ 
tially separated and orthogonal. They satisfy = 

'I'Mj, j = 1,2 where S = OyTyC is the particle-hole op¬ 
erator and C the operator for complex conjugation. The 
analytical calculation of the MBS we present in the Sup¬ 
plemental Material [38]. Note that the MBS remain at 
the zero-energy for nonzero p (see [38] for the proof). 
The probability densities of are shown in Fig. [^ for 
different values of 4>i — 4 > 2 - The result shows that their 
positions, which are defined by the locations where the 
probability density takes its maximum, are on opposite 
sides of the junction’s circumference, and vary as a func¬ 
tion of (j)i — 4>2- We denote the positions hy 9 = 9± with 
9p{9-) for 'I'Mi('kM 2 )- They are determined by the con¬ 
dition A(j) = ±7r and, thus, can be written as 


0 ± = 


(/)2 — (/)1 ± TT 
2 


(4) 


where At/) = (j)i — (j )2 + 29 is the local superconducting 
phase difference across the junction. It is consistent with 
the previous studies that Majorana fermions can be found 
in linear [Hilo] or circular [55] Josephson junction when 
the (local) difference of the superconducting phase is ± 7 r. 


(^2= 



(c) (|)i - (j) 2 = It (d) {|)i-(|)2=3;t/2 
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FIG. 2: (a)-(d) Spatial probability densities of the zero-energy 


MBS 
Eq. (|4 
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^Mif a-nd |^M 2 r, for different values of t/)i — 02 [see 


The interface of Si and S 2 is depicted by the white 
circle. The MBS are exchanged (or rotated by tt) when 
01 — 02 varies from 0 to 27r. The parameters are A = 5^ and 
fj, — 0, where ^ = hvp/Ao is the superconducting coherence 
length. 


The fact that the positions of the MBS are moved by 
changing 0i — 02 allows us to perform an adiabatic ex¬ 
change. When we change adiabatically either 0i from 0 
to 27r or 02 from 0 to —27r, the two MBS rotate by tt 
in a clockwise direction and are exchanged, as shown in 
Fig.H During this operation, they continuously evolve 
as follows [7]: 


'f'Mi —f —s'I'm2, 4'm 2 —f s'I'mi, (5) 

where s = ±1. s = 1(—1) corresponds to the change of 
01 ( 02 ) mentioned above. We note that this transforma¬ 
tion for the exchange is exact as long as the adiabaticity 
is fulfilled because the bound states are exactly at zero 
energy. 

The adiabatic change of 0i — 02 can be achieved ex¬ 
perimentally by applying a constant voltage Vj across 
the junction. For the adiabaticity, eVj should be much 
smaller than the excitation gap to the excited states. 
Then 0i—02 varies with time t as 0i—02 = 0o + 2eVjt//i, 
where 0o is constant in time. In this case, mi ^^nd 
M 2 become instantaneous zero-energy eigenstates of the 
time-dependent 'HBdG[0i(O) <(’ 2 (i)], and their evolution 
during a half rotation time Tj = T:h,/{eVj) is described 
as in Eq. § m- Equivalently, in terms of Majorana 
operators, the exchange is expressed as 


7 i(Tj) =-372(0), 72(^7) = S 7 i( 0 ). 


( 6 ) 
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Here, 7 i=i, 2 (i) = J d'^r and $(f^ = 

(d)^, —$ 1 )”^ is the four-component field operator. 

Exchange phase and transport .— We now consider a 
normal metallic tip (STM) tunnel coupled to a position 
of the Corbino Josephson junction containing two circu¬ 
lating MBS. In this situation, the coupling between the 
tip and one Majorana state becomes significant or sup¬ 
pressed as the Majorana state approaches to or leaves 
from the tip, respectively, while the other Majorana state 
is on the opposite side of the junction with an exponen¬ 
tially small coupling to the tip, leading to occurrence of 
time-dependent tunneling periodically with period Tj. In 
each period, MBS are exchanged. We will see below how 
the resulting exchange phase influences the tunneling dif¬ 
ferential conductance in this setup. 

We perform our calculation in the weak tunneling 
limit. Also, in order to focus on the low-energy regime 
where the tip-Majorana coupling dominates, we assume 
that the voltage eVj(= ttH/Tj) across the junction and 
tip voltage eV are much smaller than the excitation gap. 

Then the low-energy Hamiltonian of our setup is 

H{t) = Hn + HBdoit) + Hrit). (7) 

H]s[ = <7 ^kcl^Cka describes the normal metallic tip 

where fc, a represent momentum and spin, respectively, 
and HBdcit) = 5 / df^^r)'HBdG{t)^{t^ ■ The tunneling 
Hamiltonian between the tip and the junction is 

= d?r cl^ + H.c.]. ( 8 ) 

fc ,(7 

As we are interested in low energies, we project 
Ht onto this low-energy subspace by using ^a{r) — 
Then we obtain the time-dependent 
tunneling Hamiltonian, 


where 14 = bi,fecr(to)- Next, around t = tq, we assume 
that the magnitude of the coupling coefficient varies with 
time in an exponential manner, while the phase of the 
coupling does not change significantly; the former as¬ 
sumption is justified by the fact that the MBS show in 
Fig-i are exponentially localized in the azimuthal di¬ 
rection |38] . and the latter one is a good approximation 
provided that <C Tj, where is the tunneling du¬ 
ration. Then we have 

i7T(i) = Ee”^'‘”‘’'[^''47i(to) + H.c.]. (12) 

q,k 

The current through the tip is {I{t)) = —e{dNT{t)/dt). 
Nrit) = ^Fe tip number operator in the 

Heisenberg picture with respect to H{t). The expecta¬ 
tion value (•) is taken over a thermal ensemble of initial 
states at time t = to being in the far past with Ht = 0 . 
We switch on Ht at this time [12]. Since tunneling is 
assumed to be weak, we approximate the time evolution 
operator of H{t) up to first order in HT{t) in the inter¬ 
action picture. We then obtain the tunneling current 

= f ^k,qq'it,t') 

I fc.9.9' 

X [Gk{t,t')-Gk{t,t')]Mit,t')^, (13) 

in the limit of A“^ ^ Tj, where the operator notation 
O means an interaction picture operator. The current 
operator is given as 


Hrit) = cl^ 7j(t) -h H.C.], (9) 

where Vj^kait) is the coupling coefficient between the tip 
and 7 j(t) at time t, given by 

Vj,ka{t) = j d^r tka{r) ( 10 ) 

For further development, we consider a specific model 
for HT{t). We first define the time sequence tq = to + qTj 
of the coupling with an integer q such that if q is even 
(odd), Vi^kcritq) [V 2 ^kcr(tq)] IS maximal, corresponding to 
a situation where 71 (tq) [ 72 (^ 5 )] is located closest to the 
tip, while 72 (tq) [ 7 i(tq)] is on the opposite side of the 
junction with an exponentially small coupling to the tip. 
Moreover, since MBS are spin polarized, we can restrict 
the consideration to electrons of one spin species of the 
tip and, thus, effectively deal with spinless electrons. By 
taking nearest-neighbor couplings and writing them in 
terms of Vi^ka{to) and 71 ( 4 ) using Eqs. ([^ and ([^, we 
obtain the tunneling Hamiltonian at t = tq as 

Hritq) = HT{to) = ^[Vkcl-fi{to) + H.C.], (11) 

k 


+H.C.]. (14) 

q,k 

We note that c^(t) and 71 (t) evolve in time under Hj^ 
and HBdoit), respectively. The tunneling magnitude is 
defined as 

= (15) 

Gk{t,t') and Gk{t,t') are the tip electron Green’s func¬ 
tions |43|, and M{t,t') is the Majorana Green’s function 

M{t,t') = -471 (t) 7 i (!'))■ (16) 

This function has the information of the Majorana ex¬ 
change that is crucial in our proposal. Since the tun¬ 
neling current is exponentially small except for t = tq 
and t' = tq' due to Tk,qq'{t,t'), we can approximate 
M{t,t') = M{tq,tq') = —i{'^i{tq)^i{tq')) . In the adia¬ 
batic limit, the time evolution operator satisfying 71 (t -|- 

Tj) = UlMi)UM is m 


Um = exp 


172 ( 4 ) 71 ( 4 ) , 


(17) 
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where a = s7r/2 is the exchange phase where s = ±1 was 
defined in Eq. From the transformations 7 i(tq(')) = 

M{tq,tq>) can be obtained as 




= _ jg*’)a(9-9') 


(18) 



-2 71 -7t-a0 ttTT 27t 


eV[;i TJ‘] 


Here, 'P(to) = —*( 7 i(^o) 72 (to)) is the fermion parity of 
the junction at t = to and is equal to 1 (— 1 ) when the two 
MBS share no (one) fermion, and rj = ±1 for 7^(to) = ±1- 
By substituting Eq. (18) into Eq. (13), we can obtain the 
tunneling current (/(t)) as a function of time t and the 
tip voltage eV. 

We discuss the voltage dependence of the time- 
averaged differential conductance dl/dV [SS]. It is 
enough to consider the current for a = -k 12 and ij = 1 
as it is the same for a and 7 with opposite signs. In the 
limit of ^ Tj, we obtain 


dl 


^ sech^ 
dfcBT 


where 


and 


Zz(a) =27rr 


2kBT 


2XTj 


■ sech 




2kBT 


1 2 


X^Tj + {a - 2 ttI)^ 


Vi^{a) = ±eE — —(a — 27r/), 


h 

X' 


(19) 


( 20 ) 


( 21 ) 


with I = 0,±1,±2,.... Here, E = 2-Kp\Vk\^ is the tun¬ 
neling energy, and p is the tip density of states at the 
Fermi level. We assume the wide-band limit for the tip 
so that r becomes energy independent. The calculated 
di/dV is valid for kBT \Zo(a)\. It shows peaks at 
eE = ±{a — 2'Kl)h/Tj with height {e^/h)Zi{a)/(AkBT) 
as shown in Fig. 

Discussion and experimental feasibility .— We have 
shown that the time-averaged differential conductance 
between the metallic tip and the Corbino-Josephson junc¬ 
tion hosting two rotating MBS exhibits peaks at eE = 
±(a — 2TTl)h/Tj, regardless of the fermion parity of the 
junction. This feature results from the coherent inter¬ 
ference in the time-dependent tip-Majorana tunneling 
where the exchange phase a from the half rotation of 
MBS introduces a relative phase between the tunneling 
at time t' and t' + Tj, as shown in Eq. (18 1 . We empha¬ 


size that, different from the case of tunneling between the 
tip and a static MBS which yields a zero voltage conduc¬ 
tance peak [15], we have conductance peaks at nonzero 
voltage for zero-energy MBS due to the exchange opera¬ 
tion. The height of the conductance peaks decreases with 
increasing tip voltage. This is because when the voltage 


FIG. 3: Time-averaged differential conductance as a function 
of eE. The parameters are = 0.1 meV = 10“^?iA = 

lOfcsT = lOr. The exchange phase a = 7r/2 is used. Around 
zero voltage, the peaks emerge at eE = ±ahTy^. 


increases, the correlation time of electrons in the metal¬ 
lic tip h/ (eE) becomes much smaller than Tj, leading to 
the suppression of the interference between tunneling at 
t' and t' -I- Tj and, hence, to a decreasing peak height. 
We note that the position of the conductance peaks can 
also be understood by considering an effective Floquet 
Hamiltonian as the BdG Hamiltonian is periodic in time 
with period Tj [55] . 

We now discuss the experimental feasibility. In our 
proposal, we require thin-film superconductors whose 
Pearl length A is much larger than the values of R and 
f. Thin-film NbN superconductors have A ~ 30 pm, 
^ ~ 3.5 nm, and superconducting gap Asc ~ 3 meV [44] . 
Assuming the proximity-induced gap Aq ~ 1 meV, the 
excitation gap of our Corbino-Josephson junction can be 
estimated by Eg ~ ^j2/S.QhvF/R ~ 0.63 meV for R = 
5Hvf/Ao and p = 0, where vf = vfAI/{p'^ Aq) [5]. 
Next, in order to have adiabatic rotation of the MBS and 
a short tunneling duration, Tj should be long so that 
h/Tj{= eVjj'K) <C mm{hX,Eg}, but sufficiently short 
so that Tj <C tqp where tgp is the quasiparticle poi¬ 
soning time that is in the range of 0 . 1—1 ps [IS]. Esti¬ 
mating mm{hX, Eg} ~ 0.63 meV and h/tqp ^ 6.6 neV, 
gives h/Tj ~ 10^ peV for adiabatic interference with¬ 
out quasiparticle poisoning. In order to observe clearly 
the conductance peaks shown in Fig. [^ the temperature 
should be smaller than h/Tj which provides an experi¬ 
mentally feasible temperature range of 10—10^ mK. We 
finally note that the MBS are robust against weak disor¬ 
der; i.e., they remain at zero energy [55], and, hence, the 
disorder has no significant effect on our results in Fig. [^ 

In conclusion, we have proposed an experimental setup 
where the exchange phase of mobile MBS can be probed 
as a result of their interference effect. Different from the 
preceding detection schemes relying on charge neutrality 
of zero-energy MBS, our setup allows us to identify the 
nontrivial MBS statistics, providing alternative routes of 
MBS detection. 
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Supplemental material for “Detecting the Exchange Phase of Majorana Bound States 
in a Corbino Geometry Topological Josephson Junction” 

Sunghun Park^ and Patrik Recher^’^ 

^ Institute for Mathematical Physics, TU Braunschweig, 38106 Braunschweig, Germany 
^ Laboratory for Emerging Nanometrology Braunschweig, D-38106 Braunschweig, Germany 

In this Supplementary Material, we analytically derive Majorana zero-energy states for /r = 0 and prove that they 
are still at zero energy for finite fx. Next, we provide the details of the calculation of the time-averaged tunneling 
current I in the main text, and an alternative explanation for the result shown in Fig. 3 using Floquet theory. Finally, 
we discuss disorder effects. 


A. Derivation of Majorana zero energy states 


In this section, we first calculate Majorana zero-energy states, 'I'mi and 4 'm 2, for ^ = 0 by solving the BdG equation 
for the Corbino geometry. From Eqs. (1) and (2) in the main text, the BdG equation for /i = 0 is given by 


/ vpd ■ p 

[A*{r,9)l 




( 22 ) 


This decouples into two sets of equations, one for spin up and one for spin down. We consider only the equations 
with spin down as those for spin up give diverging solutions which are not physical. Due to the absence of rotation 
symmetry of the system, the wave function in each region should be given in the form of a superposition of different 
angular momentum eigenstates. Hereafter we will measure energy and length in units of gap magnitude Ag and 
superconducting coherence length ^ = Hvf/Nq, respectively. With these units, the wave function for the region of 
0 < r < i? is expressed as 


OO 


Here, are coefficients to be determined below, and 0) is defined as 




( 0 

V 0 


\ 

/ 


(23) 


(24) 


where Jmiir) is the Bessel function of the first kind that is regular at the origin r = 0. For r > R, the wave function 
is found as 


OO 

n= —OO 


(25) 


with coefficients 5„, and 


0 


^nir,9) = 




ir) 

0 


(26) 


with Hn \ir) the Hankel function of the first kind that goes to zero as r —)■ oo. In order to obtain solutions, we 
determine the coefficients and by matching the wave functions 4'<(r, 0) and 'I'>(r, 0) at r = R, leading to 


ai+ie-*^i/2jz+i(zR) = RRi{iR) 


(27) 
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and the following recurrence relations can be derived, 


{^R)Ji+i {^R) 
6 ,,^, = Hi,{iR)Ji,+2{iR) 


ai, 


bi', 


( 28 ) 


where I and I' range over all integers. Eqs. (271 and (28) allow us to express the coefficients and &/'=odd («/=odd 

and &z'=even) In terms of gq (a_i). This fact indicates that we can construct two independent solutions, denoted by 
'I'l and tI/ 2 , 


^'i=0(-R-O + 0(r - i?) 

m—even n—odd 

vl/2 = 0(i?-r) ^ am'^<+eir-R) ^ bn^>, (29) 

m=odd ra=even 


with one unknown coefficient ag (u-i) for Ti ('^ 2)5 where 0(x) is the unit step function. According to the particle- 
hole symmetry, if 'ki is a solution, then STi is also a solution, where S = <^yTyC is the particle-hole operator 
with C the operator of complex conjugation. This imposes a constraint '^2 = S'?!, and hence a_i = —Uq. We 
choose oq = and a_i = —where N is the normalization constant. The following 

combinations of 'I'l and ^['2 give two MBS, 


1 % 

^ Ml = -^{'^1 + ^2), ^ M 2 = ^{'^1-^2), 


(30) 


which are plotted in Fig. 2 in the main text. 

We show that 'I'mi and centered at 0 = 0±, respectively, are exponentially localized in the azimuthal direction 

to justify the assumption in Eq. (12) in the main text that the coupling between the tip and MBS is exponentially 
small except for t = tq where tq = tg + qTj with an integer q. We first solve an effective low-energy Hamiltonian which 
is the k ■ p Hamiltonian obtained by Fu and Kane in Ref. m incorporating the spatial variation of a superconducting 
phase difference due to two flux quanta in the Josephson junction, 


'Heff = 


fAgCOS 



(31) 


where v = vf[cos{p,W /hvp) + {Ag/fj,)sm{pW/hvF)]AQ/+ Aq), p is the chemical potential, W is the distance 
between two superconductors, and Ai^ = — (^2 + 20 is the phase difference of the Josephson junction. The effective 

Hamiltonian, whose bases are two branches of counter-propagating helical Majorana states after integrating out the 
transverse (radial) degree of freedom, is a good approximation in the limit of i? where = Hv/Ag. To calculate 
approximate zero energy solutions localized at 0 = 0p with p G —}, we linearize cos(A(/)/ 2 ) around A(j) = pii, 


cos 



= cos 


^pTT -I- 2d0p^ 


= —p sin((50p) 
« -p SOp 


(32) 


for small 69 p where d 9 p = 9 — 9 p. Then we obtain the exactly solvable Hamiltonian 


/ -ifde -ipAg 69 p\ 

\ ip AgS 9 p ib§de j ’ 


(33) 


thus allowing us to calculate the zero energy solutions or approximate Majorana states which are exponentially 
localized in the azimuthal direction, given by 


^approx ( 0 ) 


R 


R 


(0 - 9pf 



(34) 









1.5 


in 

C 
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FIG. 4: Probability densities of MBS as a function of azimuthal angle Q for the values of d± = ±7r/2, R = 5^' and /r = IF = 0. 
Probability densities of in EQ- (361 nre drawn by the solid line, and those of approximate solutions in 

Eq. (351 are drawn by the dashed line. 


and probability densities 


I ^approx (0)12 



(35) 


To compare these approximate solutions with the exact MBS given in Eq. (30), we integrate out the radial degrees of 
freedom of and '^M 2 {r, 0 ), 


l'I'MifM2)Wl" = / drr\^MiiM2){r,e)f. (36) 

Fig. 1^ shows the comparison between 'I'®='mi,M 2 (^) iE'pL+™(^)j exhibiting a quantitatively good agreement 

between them. This leads us to the conclusion that MBS found in Eq. ( [M| ) are exponentially localized in the azimuthal 
direction, and hence the form of the tunneling Hamiltonian of Eq. (12) in the main text is a good approximation. 


B. Zero-energy Majorana bound states for non-zero chemical potential 

We prove that Majorana bound states remain at zero energy for finite fj, by using a Taylor expansion around fi = 0. 
We start the proof by rewriting Eq. (1) in the main text as 

(37) 

where are the two lowest eigenvalues which are assumed to be differentiable functions with respect to /r, and 

'^j=i, 2 {R) are corresponding eigenstates. In the limit of /r — >• 0, Ej{fj,) = 0 and ^'j(/i) become two Majorana bound 
states, 'I'i(O) = and 'I'2(0) = 5 'm 2- Throughout this section, we use the notation I'll) = Sldl) as the particle-hole 
transformed counterpart of I'll) where S = TyayC is the particle-hole operator. We note that ('F'l'E') = ('F'l'E')* by the 
anti-unitary property of S. And I'll) = I'll) if I'll) is a Majorana state. Below, we will show that Ej{^) = 0 also for 
finite /i. 

The Taylor expansion of the BdG equation in terms of /r leads to 

’HBdG(M) =’^BdG(O) —/rr^, (38) 

Ej{^) = Ej^o + Ej^i^ + + ..., (39) 

^j(m) = + • ■ •, (40) 








9 


where 


d-Eifi) 


fi—0 




(41) 


^—0 


and N is the normalization constant. We note that Ej^o = 0 and 'I'jo = '^Mj- By substituting these series into 


Eq. (371 and equating the terms in each power of /i, we obtain equations for zeroth-order, first-order, and so on, each 


of which is discussed below. 

The equation for the zeroth order of /i is given by 


’HBdG(0)4'j_o = = 0, 


(42) 


which is trivial. 

Next, we consider the first-order equation given by 


If we multiply by (Tj^ol, it yields 


’HBdG(0)4'j-,i — T’z'I'yo — 


(43) 


(«',.o|HBdG(0)|4',.i) - (vI^,.o|r.|4-,.o) = E,, 1(4-,,o|4/,,o). (44) 

The first term on the left side is zero by virtue of Eq. ( [4^ . Then the second term is real, and we have 

(«'i,o|r.|4'„o) =(4-,-o|r,|4-,,o)* = ($,, olSr.lT,, o) 

= - (^..olr.l^-j^o) 

=0, (45) 


resulting in Ej^i = 0. Here we used anti-unitarity of and anticommutation relation 
Ej^i = 0 leads to 


0. Eq. (43) with 


If we apply 5 to Eq. (46) 


allowing us to write 


'HBdG(0)'I'j.l = 

('I',',o|r.|4-,.o) =0 Vj,/g{1,2}. 


'HBdG(0)4',,i =t^4',-,o, 

= 'HBdG(0)'I-,-,i, 




(46) 

(47) 


(48) 


(49) 


which will be useful to prove Ej 2 = 0 below, where a,,i and /3,,i are constants. 

We can do the similar procedure as above to get i ?,,2 = 0. The /r^-order equation with £’,,o = Ej^i = 0 is given by 

^'HBdG(0)4',-,2 - (50) 

Multiplying (4'j,o| to this equation gives 

-('I^,.o|t.|4-,-i)=E,-2(4-,.o|%o). (51) 


Using Eqs. (47) and (49), 


(«'j,o|t.|4',,i) = (4',,o|t.|4-,,i)* = -(4',,o|t.|4-,,i) 

= - aj,i(4'j,o|Tz|4'i,o) - /3i,i('I'i,o|'rz|4'2,o) 

= 0 , 


(52) 











10 


and we get Ej ^2 = 0. From Eq. (50), 


■^'HEdoiO)'^ j,2 = Tz'SjA 

('f'/.o|r.|vI/,-1) = 0 Vj,/g{1,2}, 

'Sj,2 = ^'i,2 + + Q;j,2'f'i,o + /3j, 2 ^ 2 , 0 , 


(53) 

(54) 

(55) 


where aj ^2 and /3j_2 are constants. 

We now verify that Ej^n = 0 for all positive integer n by induction on n. It is assumed that the following conditions 
from the ^"-order equation hold: 


Ej^k = 0 for 0 < k < n, (56) 

'HBdGiO)'^j,n = riTz'^j^n-l (57) 

(4'i'.o|Tz|4'j,fc') = 0 Vj, j'G {1,2}, and for 0 < A:'< n - 1, (58) 

n—1 

'^j,n = '^j,nE Cj! ^k'^ j' ,k', (59) 

j' = l,2 k'=0 


where Cji^k' are constants. Since we have shown above that these conditions are fulfilled for n G {1, 2}, it is enough 
to show that if they hold for n = uq, then also n = hq + 1 holds. Using Ej^k = 0 for 0 < fc < no, the -order 
equation is given by 


(no + i)!^BdG(0)4/„„„+i ^ F;,,„„+i'h,.o. 


If we multiply (4'^ o| to this equation, we have 


- A(^j-.o|T.|4',,no) = 7- ^F;,-„„+i(4',-o|4',.o). 


no! 


(no -f 1)! 


Using Eqs. (57) and (58) with n = no, ('I'yo|T^|'l'j,no) i’a the left side can be calculated as 


— ('I'yO I'T'z |) 
"0-1 

j' = l ,2 fe '=0 

= 0 , 


and Ej^no+i = 0- From Eq. (60) we get 


77BdG(0)'I'i."o + l = ( 1*0 + l)l‘z'f'i,"o> 
('l'/.o|r.|vli,,„„)=0 Vj,/g{1,2}, 

"0 

^t,»o+i=^i,»o+i + y~! ^j',k'^j',k'■ 

j' — l,2 k'—O 


Therefore, if Eqs. (56l-(59) hold for n = no, then they also hold for n = no -I- 1. 


(60) 


(61) 


(62) 

(63) 

(64) 

(65) 


C. Time-averaged tunneling current 


In this section, we provide the details of the calculation of the time-averaged tunneling current I in the main text. 
It is obtained from the average of {I{t)) over a half rotation period. 


I 


1 



dt{I{t)), 


( 66 ) 
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in the limit of very large integer Q so that tg = to + QTj ^ to. For a time t G [tq — Tj/2,tQ +Tj/2], (/ (t)) is roughly 
proportional to the tip-Majorana tunneling magnitude {I(t)) oc Thus {I{t)) is maximum when t = tg, and 

is exponentially small of the order of for t = tg ± Tj/2. 

We first calculate the tunneling current in the time interval [tg — Tj/2,tg + Tj/2]. From Eq. (14) in the 

main text, we obtain 


m) 


2e |VfcpA 

A2 + (efc + eE)2/fi2 


[-(1 - nF{ek))A{t, a) + nF{ek)A{t, -a)], 


(67) 


where 


A(t,a) = e^At-tQ) 

for t e [tg - rj/2,tg], and 


■g-i 

) ^ina^-i{ek+eV){t-tQ-„)/fi 


+ 




2A 


+ c.c. 


( 68 ) 


A{t,a) = -e-^At-tQ) 


■g-i 

.71 = 0 


^^±A^±±ey)/^^^ga^-i{ek+eV){t-to)/h 

2A 


(69) 

for t G [tg,tg +Tj/2]. A{t,—a) is obtained by replacing a by —a from A{t,a). In this calculation, we neglected 
the terms of the order of that give an exponentially small contribution to A{t,a) in our limit of the short 

duration of the tunneling A“^ AiTj. Then I can be expressed as 


I 


2e |VfcpA 

/i2 ^ A2 + (efe + ey)2/;i2 


f 

\ rtQ + A~ 


1 ft'Q+A 

< -(1 - UFiSk)) 

— dtA{t,a) 

[Tj J 

+ UF^Sk) 

— dtA{t,-a) 

Aq-A j 


To calculate / we have to calculate the time averages of A{t, a) and A{t, —a). It is given by 


(70) 


\ r*Q+'^ 

— dtA{t,a) 


1 


2A 

A2 + (efc + eE)2/fi2 


Q 

^-K^k+eV)nTj/h 

n——Q 


-A + l{£k + eV)/h ,ag i(^^+^v)gTj/n 

A2 + (efc + eE)2/fi2 


+ c.c. 


(71) 


In the limit that Q goes to infinity, the summation term over n on the right side can be written as 


Q 


lim 

Q—>-oo 


'y ^ ^ina^-i{€k+eV)nTj/h 
n——Q 


=2f J2 ^ 


/ = —OO 


{ek + eV)Tj 
a -r-27ri 


2Fh 




—Sk — eV + 


(a — 2fI)H 


(72) 


where I = 0, ±1,±2,..., while the remaining terms that possess factors can be neglected because of their 

rapid oscillation as a function of k for large Q 


E 


^±iekQTj/h 


k 


0 . 


(73) 
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Thus the time averages of A(t, a) and A(t, —a) (obtained by replacing a by —a in A(t, a)) have the forms 

(a — 2TTl)h 


1 /■‘Q+^ j 4, X 2A 2TTh ^ ^ 

X^ + ie^ + eVy/h^Tf ^ 




-£k - eV + 


Tj 


1 /■*« + ■ 


dtA{t, —a) = 


2A 


2TTh 


A2 + (£fc + evy/h^ Tj 


E ^ 


r=-c 


2A 


2Trh 


A2 + {Sk + evy/h^ Tj 


E ^ 


-£k -eV - 

-£k - ey - 


(a + 2nl')h 
(a — 2TTl')h 


(74) 


In the last line of the second equation, we changed I' by —V for convenience. If we put these into Eq. (70), it follows 
that 




2A 


A2 + (efc + eE) V;i2 


- (l-ni.(£fe)) S 


l = — oo 


—Sk — eV + 


(a — 2Trl)ti 


+ nFi£k) Y ^ 


l' — — C 


-£k -eV - 


Tj 

{a — 2'Kl')h 


Tj 


By replacing the discrete sum over k by an integral with the density of states of the metallic tip p{£k) 

Y^ tiefep(£fc). 


(75) 


(76) 


Assuming a wide-band limit with an energy independent tunneling rate T = 27rp(efc)|Vfcp, we obtain 

^ ZAJ- ] "* 

. N 97rr - - - 

h 


— 6 V—> 

7 =jy: 2 ,r 


X^Tj + {a-2Triy_ 


[—npieV — {a — 2TTl)ti/Tj) + np^—eV — {a — 2'Kl)h/Tj)] 


= T Y [-^F{Vy{a)) + npiy (a))] , 


(77) 


where Zi{a) and Vj^^a) are defined in the main text. It is now straightforward to evaluate the time-averaged 
differential conductance dl/dV. Since only the Fermi-Dirac distributions have the voltage dependence, it is enough 
to compute 


dnF{Vi^{a)) d 


1 


dV dE 1-I-exp [V)^(a)/(fcBr)] 


= T 




sech^ [E±(a)/(2fcBr)] . 


Therefore, we get 


di 


ai e 
PV ~ ~h 2^ 


dV h AksT 


V)+(«) 

2kBT 


■ sech 


vn^) 

2kBT 


(78) 


(79) 


D. Floquet description of the two rotating Majorana bound states 

In the previous section, we showed that the time-averaged differential conductance has peaks at eV = ±(a — 
2'Kl)h/Tj. Here, we provide an alternative derivation of the result by considering an effective Floquet Hamiltonian 
for the two rotating Majorana bound states in our Corbino-geometry Josephson junction. The Floquet description 
is valid as the time-dependent BdG Hamiltonian in Eq.(7) in the main text is periodic in time with period Tj: 
i7BdG[A(/)(F)] = 77BdG[A<(i(t'-I- Tj)] where IXcj) is the phase difference across the junction. The time evolution operator 
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of the Hamiltonian from t' to t' + Tj in our low-energy subspace is Um in Eq- (17) in the main text, and the effective 
Floquet Hamiltonian Hp |S2) is obtained by using 

e-^HFTj/h ^ (gQ) 


which yields 


Hp = i^72{to)ji{to) + (81) 

where e = ah/Tj and integer 1. The first term is reminiscent to the Hamiltonian for two coupled static MBS with 
hybridisation e (which shows conductance peaks at eV = ±e), and the second term shifts energy levels. As a result, 
tunneling into an effective system described by Hp gives conductance peaks at eV = ±e + 2TTlH/Tj which is producing 
our features above. 


E. Weak disorder effect 

In this section, we show that two Majorana bound states in a Corbino-geometry Josephson junction are robust 
against weak disorder on the surface of a topological insulator if we stay in the subspace spanned by the two Majorana 
bound states, i.e., they remain at zero energy. It can be shown by using a Taylor expansion in the same way as in 
Section B. The surface of a topological insulator with a disorder potential is described by 


H = vpa • p + V{f), 


(82) 


where V (r) is the disorder potential. Then the BdG equation is given by 


/ vpd ■ p A{r,6)I \ . f XV{f) 0 \ 

l^A*(r,6»)I -VFa-p)^\ 0 J 


mr,e) 




(83) 


Here a perturbation parameter A is introduced to keep track of the order of the disorder potential in a Taylor expansion, 
and 'I'j=i_ 2 (’"j E' are eigenfunctions and eigenvalues, respectively, which are assumed to be differentiable 

functions with respect to A. Like in the case of the calculation in Section B, we do a Taylor expansion for 2 (''’) ^) 
and Ej in A, 


^=^.o + AE',i 


4'' = 5'' n 

3 3,0 


A2 


3! 


A2 

^ j.i ^ 2! 3! 


(84) 

(85) 


Now we substitute these expansions into the BdG equation and equate terms of the same order in A. All calculation 
details are the same as in Section B, except that p and —Tz in Section B are replaced by A and V(r)Tz where Tz is 
the Pauli matrix in particle-hole space, and we get 


E'.fc = 0, (86) 

for all j G {1, 2} and nonnegative integer k. Therefore we conclude that two Majorana bound states are robust against 
a weak disorder potential, and hence the results shown in Fig. 3 in the main text can be detected even in the presence 
of weak disorder. 
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